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Abstract. The first measurement of the dielectron continuum at RHIC energies was performed by the 
PHENIX experiment for Au+Au collisions at y / sjvjv= 200 GeV. Mass spectra for different centralities are 
presented and compared with the expectations from hadron decays. 



1 Introduction 

Electromagnetic probes are ideally suited to investigate 
hot and dense matter produced in high energy heavy ion 
collisions because they do not undergo strong interac- 
tions and thus probe the time evolution of the collision. 
The dielectron continuum is rich in physics. In the low 
mass region Dalitz decays of light hadrons and direct de- 
cays of vector mesons, which might be modified in the 
medium contribute to the spectrum. A number of ex- 
periments (E325[l], DLS0, CERES and more recently 
NA60 4 ), independent of bombarding energy, have ob- 
served an excess of dielectron yield over the hadronic sources 
by a factor 2-3 for masses between 0.2 and 0.8 GeV/c 2 , 
when going from proton to heavy ion induced reactions. 
This enhancement has been interpreted as thermal radia- 
tion from pion annihilation in the hot fireball largely me- 
diated by light vector mesons p, oj and <p. Among these, 
the most important contribution arises from the p meson, 
due to its strong coupling to the 7T7t channel and its short 
lifetime of only 1.3 fm/c which, in contrast to the longer 
living lu (23.4 fm/c) and </> (44.4 fm/c), makes it very suit- 
able to probe in-medium modification close to the QCD 
phase boundary. 

In the intermediate mass region between the cj> and the 
J /\P vector meson the dominant contribution arises from 
correlated charm production. This region has been pro- 
posed as an interesting candidate to search for thermal 
radiation emitted as a blackbody radiation from the QGP, 
since its contribution could be comparable to that of the 
charm and would be dominant to higher pt Earlier 
measurements |H] indicate an excess of yield in nuclear 
reactions with respect to the rate expected from elemen- 
tary reactions. More recently NA60 |7] has disentangled 
the prompt dimuons from the pairs originated from open 
charm semi-leptonic decays and demonstrated the prompt 
origin of the excess yield. 

Although correlated e + e~ pairs are rare, the 0.24 nb _1 



collected by PHENIX for Au+Au collisions at ^snn= 
200 GeV in 2004 provide a significant sample to investi- 
gate the dilepton continuum. In the following, the analysis 
of 800M events is presented. 



2 Di-Electron Analysis 

The PHENIX detector design is optimized for electron 
measurements. It combines an excellent mass resolution 
(1 %) with a powerful particle identification achieved by 
matching the reconstructed tracks with the information 
from a Ring Imaging Cherenkov detector (RICH) and an 
Electromagnetic Calorimeter (EMC). Electrons are iden- 
tified by requiring at least 3 phototubes matched in the 
RICH and by correlating the energy measured in the EMC 
and the momentum, parametrized in terms of E ~ p . Pair 
cuts are applied to avoid sharing of detector hits, tracks 
which are parallel in the RICH are rejected by a cut on 
the angular difference. Whenever encountering such a pair, 
the complete event is rejected. 

Pairs are created by combining all electrons with all positrons 
in one event. The overwhelming yield of these pairs is un- 
physical. A statistical procedure is used to determine this 
combinatorial background. Since the PHENIX acceptance 
is different for like and unlike sign pairs, the combinato- 
rial background is computed with a mixed-event technique 
by pairing un-like sign tracks from different events within 
similar topology (i.e. vertex position and collision central- 
ity). The like-sign pairs prove that the shape of the back- 
ground is reproduced with high precision by the mixed 
event technique, since the ratio of the like-sign spectra for 
real and mixed events deviates from 1 by less than 0.1 %. 

Four different methods have been tested to normalize 
the background distributions. One normalizes the number 
of mixed events to the number of physical events. The sec- 
ond one relies on the mean number of single electron tracks 
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which contribute to uncorrelated pairs (iV_| ) = rTjl • nZ. 

Both these methods normalize directly or indirectly to the 
product of the average electron and positron multiplicities, 
which is equivalent to the assumption that the sources of 
electron and positron tracks are independent. Unless the 
probability density is a Poisson distribution an additional 
correction factor 



and are given by 
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is introduced to account for the primary multiplicity dis- 
tribution. Moreover in heavy ion collisions the probability 
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Fig. 1. Unlike sign mass spectrum: foreground (black), back- 
ground (red), subtracted spectrum (blue), with systematic er- 
ror marked in green (all the systematics except for the back- 
ground) and yellow (systematics of the combinatorial back- 
ground). Colors in the online version. 



to produce N pairs depends on other variables which will 
cause additional correlations in the number of pairs pro- 
duced in an event. The most relevant are multiplicity cor- 
relations due to the finite width of centrality classes, which 
are characterized by the number of participants N part . The 
corrections have been evaluated with a Monte Carlo based 
on the Glauber model parameterization of N part and N co u 



(N+)(N-) 



(2) 



Similar studies have demonstrated that correlations due 
to the finite size of vertex classes are negligible, consider- 
ing the low multiplicity of the electron sample. 
However, due to the fact that electrons and positrons are 
always created in pairs, the unlike sign background is the 
geometric mean of the like sign backgrounds, indepen- 
dent of the primary multiplicity distribution (N^ ) = 

2 • ^N ++ N This normalization turns out to be statis- 
tically correct and insensitive to multiplicity correlations 
like the first two. However, due to the presence of a small 
signal in the like sign yield due to double conversion and 
double Dalitz decays, the geometrical mean of the like sign 
yields overestimates the combinatorial background. This 
problem can be repaired by normalizing the background 
to the foreground in the like-sign distributions and exclud- 
ing the low-mass region where a signal is present; since the 
mixed events produce the proper rate of like and unlike 
sign pairs, the same factor is then applied to the unlike 
sign distribution. 

Any cut applied on a pair can potentially distort the rela- 
tive rate of like and unlike sign yield. The effect of pair cuts 
on the background normalization has been studied with a 
Monte Carlo simulation and survival probabilities for like 
and unlike sign pairs have been computed as a function of 
the event topology. Therefore an absolute (k±) or relative 
(k = -^==== =) pair cut bias correction has to be applied 

in the first two or in the last two cases respectively. 
Finally all the normalizations agree within 0.5 %. The em- 
pirical normalization of the like-sign pairs was chosen for 
the final results and a systematic uncertainty of 0.25% 
was assigned. 

After removal of the combinatorial background, the con- 
tribution of photon conversion was removed, cutting on 
the orientation angle of the pair in the magnetic field. Fi- 
nally the spectra are corrected for efficiency such that the 
data represent the dielectron yield produced in a collision 
for which the electron and the positron are both in the 
detector acceptance. 



3 Results 

Figure ^ shows the invariant mass spectrum of the mea- 
sured dielectron pairs, the background, and the subtracted 
yield with uncertainty: the error bars in the upper panel 
correspond to all the systematic errors, which are domi- 
nated by the background normalization. The lower panel 
shows separately the uncertainty arising from electron iden- 
tification, efficiency, acceptance and run-by-run fluctua- 
tions, estimated around 22% represented by the point-to- 
point green band, and the very small uncertainty on the 
combinatorial background (0.25 %) represented by the yel- 
low band around zero, which becomes the dominant source 
of the signal uncertainty due to the low signal to back- 
ground ratio. 
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Figure |21 shows the data with the total systematic error 
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Fig. 2. Data compared to a cocktail from all the hadronic 
sources. The systematic uncertainty of the data is shown in 
yellow. Colors in the online version. 



A significant signal is left over the full mass range, cor- 
responding to an integral of 1.8 • 10 5 pairs, out of those 
15,000 above the ir° mass. The data have been compared 
to a cocktail of hadron decay sources. The pion spec- 
trum used as input is determined by a parametrization 
of PHENIX charged and neutral pions. The spectra of the 
other mesons are determined from the pion spectrum by 
rriT scaling i.e. using the same Hagedorn parametriza- 
tion as for pions and replacing px with 
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The normalization of the yield relative to tt° is determined 
by the asymptotic ratio at high px (5 GeV/c is used here) 

ei cm 

The systematic error, depending on the pion yield and 
the relative cross section of the other contributions, varies 
from 10% to 25%. 

An additional source of dielectron pairs, which becomes 
the dominant continuum contribution for invariant masses 
above 0.5 GeV/c 2 is the correlated charm production. This 
contribution has been simulated with PYTHIA, scaling 
the p+p equivalent cc cross section of 622±57±160 /ibarn 
to the number of minimum bias Au+Au binary collisions 
9 . However it is worth to recall that the single electron 
measurement shows a suppression at high pt which 
increases with centrality. The implications for a dielectron 
invariant mass are not straightforward, since the invariant 
mass implies knowledge both of the momentum of the two 
electrons and the opening angle of the pair. The charm 
curve here is meant more as an indication than as a real 
quantity to compare with. The systematic error on the 
charm is therefore meant as an error on the cc cross sec- 
tion only, not on the shape. The calculated electron pairs 
from charm as well as from the cocktail, have been filtered 
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Fig. 3. Data (systematic uncertainty in yellow) compared to the 
cocktail (systematic uncertainty in cian) and theoretical predic- 
tions, where a rho spectral function is introduced, in vacuum 
(red), with a dropping (blue) or melting (green) scenario. Col- 
ors in the online version. 



compared to the cocktail. The data show a good agree- 
ment with the cocktail over the full mass region. The u> 
-and (f> resonances are not fully reproduced, most likely be- 
o cause of the combined effect of mass resolution and low 
signal-to-background ratio. The data overshoots the cock- 
tail in the region 0.3-0.8 GeV/c 2 . The systematic uncer- 
tainty however does not allow any strong conlusive state- 
ment. 

In FigureOlthe data have also been compared to the theo- 
retical predictions |12|.[5]. where the e + e~ invariant mass 
spectrum dN ee /dM has been calculated using different 
in-medium p spectral function and an expanding ther- 
mal fireball model. Although the systematical uncertainty 
does not allow any claim of a significant deviation from 
the known expected sources, it is interesting to observe 
that the data are even above the predictions that include 
in-medium modifications in the same mass region where 
other experiments quantified the same effect. 



4 Centrality dependency 

Despite the low statistical significance of the dielectron 
signal, the centrality dependency of the dielectron contin- 
uum has been studied. The centrality of the collision is 
determined using the combined information of the Beam- 
Beam Counter Detectors and the Zero-Degree Calorimeter 



4 



Alberica Toia for the PHENIX Collaboration: Title Suppressed Due to Excessive Length 



0-20% Au+Au @ = 200 GeV 20-92% Au+Au @ \Js = 200 GeV 




Fig. 4. Data compared to a cocktail from all the hadronic sources and the charm contribution and ratio for centrality classes 
0-20 % (right) and 20-92 % (left). The systematic uncertainty of the data is shown in yellow. Colors in the online version. 



on particle multiplicity, and translating this information 
into a centrality notion with the Glauber model. The event 
mixing as well as the normalization of the combinatorial 
background has been performed in each centrality class 
separately, always using the like-sign spectra as a cross 
check with an agreement better than 0.1 %. Cocktails from 
hadronic sources have been prepared from the input pion 
distributions in each centrality class [5]. 
Figure 0] shows data from two centrality classes: central 
(0-20%) and semiperipheral (20-92%) compared to the 
relative cocktails, absolutely normalized. The more cen- 
tral collision show a significant enhancement in the mass 
region between 0.3-0.8 GeV/c 2 which is not present in the 
more peripheral event sample. 

The centrality dependency is also studied via the ratio of 
yield in different mass regions with respect to the ir° yield. 
The production of pions in fact scales with the number of 
participants, with only a slight increase with centrality. If 
an in-medium enhancement of the dielectron continuum 
yield exists, it could in first order arise from 7T7t or qq 
annihilation. If so the yield should increase faster than 
proportional to the number of participants N part . Figure 
|S] shows the ratio of different mass regions to the pion 
Dalitz region (0-100 MeV/c 2 ) as a function of the number 
of participants. The left hand plot shows the region be- 
tween 150-450 MeV/c 2 where the increase as a function of 
number of participants is faster than linear. In contrast, 
charm production is known to follow binary scaling [S]. 
Therefore the yield in the charm region compared to the 
pion region should increase as the ratio ^ c °" ■ The right 
hand plot of figure El shows the region between 1.1-2.9 
GeV/c 2 , where the rather constant ratio does not reflect 
the expected behavior. 



5 Conclusions 

Thanks to the high statistics collected in 2004, the first 
measurement of the dielectron continuum at RHIC en- 
ergy has been performed by the PHENIX experiment. 
The dielectron mass spectrum is higher than the expecta- 
tions from hadronic sources and shows better agreement 
with theoretical predictions that include modified spec- 
tral function. Sampling the data in different centrality 
classes allows to isolate the excess of dielectrons in the 
most central collisions, while more peripheral collisions 
show a good agreement with the hadronic cocktail. A com- 
parison of the continuum yield in the mass region 150-450 
MeV/c 2 to the pion yield shows an interesting dependency 
of the enhancement as a function of the number of partic- 
ipants. The signal-to-background ratio however limits the 
statistical significance of these results. PHENIX is cur- 
rently developing a new hadron blind detector system to 
address this issue 1 13). 
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Fig. 5. Ratio of integral of the mass region 150-450 MeV/c 2 (left) and 1.1-2.9 GeV/c 2 (right) with respect to the ttO yield (0-100 
MeV/c 2 ). The systematic error, which depends on the combinatorial background in the region 1.1-2.9 GeV/c 2 is indicated by 
the yellow band. Note the different scale in the two plots. Colors in the online version. 
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